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Abstract 

We consider tlie effect of interactions on electrons confined to two dimen- 
sions at Landau level filling v = 2, with the specific aim to determine the 
range of parameters where the fully polarized state is stable. We calculate 
the charge and the spin density collective modes in random phase approx- 
imation (RPA) including vertex corrections (also known as time dependent 
Hartree Fock), and treating the Landau level mixing accurately within the 
subspace of a single particle hole pair. It is found that the spin wave exci- 
tation mode of the fully polarized state has a roton minimum which deepens 
as a result of the interaction induced Landau level mixing, and the energy 
of the roton vanishes at a critical Zeeman energy signaling an instability of 
the fully polarized state at still lower Zeeman energies. The feasibility of the 
experimental observation of the roton minimum in the spin wave mode and its 
softening will be discussed. The spin and charge density collective modes of 
the unpolarized state are also considered, and a phase diagram for the u = 2 
state as a function of rg and the Zeeman energy is obtained. 
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I. INTRODUCTION 



The interplay between the electron's spin degree of freedom and the inter-electron in- 
teraction has been of interest in the condensed matter physics, in particular in the jellium 
model where the neutralizing background is taken to be rigid and uniform. The relative 
strength of the interaction is conventionally measured through the parameter rs which is 
the interparticle distance measured in units of the Bohr radius as = eh^ /me^ , e being the di- 
electric constant and m the band mass of electron. For electrons confined to two dimensions, 
which will be our focus in this paper, we have 

rs^^^, (1) 

as 

where p is the two-dimensional density of electrons. The interaction strength is enhanced 
relative to the kinetic energy as the system becomes more dilute, i.e., as rs increases. It has 
been predicted that as rs is increased, the electron liquid eventually becomes spontaneously 
spin polarized to gain in the exchange energy, ultimately going into a Wigner Crystal at 
rs ~ 37 i. 

The two-dimensional electron systems, obtained experimentally at the interface of two 
semiconductors, constitute an almost ideal realization of the jellium model for several rea- 
sons. Samples with mobility in access of 10 million cm^/Vs are available f^, minizing the 
effect of disorder. Furthermore, the density of electrons, which controls the strength of the 
interaction relative to the kinetic energy, can be varied by a factor of 20 P,^. We will con- 
sider electrons in the presence of a magnetic field, specifically, at filling factor z/ = 2, which 
is a particularly clean test case for the kind of physics in which we are interested. There are 
three relevant energy scales here: huc = heB/mc is the cyclotron energy, Vq = e^/e/o is the 
typical Coulomb energy, Iq = ^hc/eB being the magnetic length, and Ez is the Zeeman 
splitting energy which is the Zeeman energy cost necessary for a single spin flip. (To an 
extent, the Zeeman and the cyclotron energies can be varied independently by application 
of the magnetic field at an angle; while the former depends on the total magnetic field, the 
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latter is determined by the normal component only.) The ground state is known in two 
limits. When hujc dominates, the ground state is a spin singlet, with 0| and 0| Landau 
levels fully occupied. On the other hand, when Ez is the largest energy, the ground state 
is fully polarized; when the Coulomb interaction is not strong enough to cause substantial 
Landau level mixing (i.e., in the limit of Ez » huc >> Vc), the ground state has Oj and 
If Landau levels occupied. When the ground state is described in terms of filled Landau 
levels, we will denote it by (A^t- ^i); where A^t and Nl are the numbers of occupied Landau 
levels for up and down spin electrons. The possible filled Landau level states at z/ = 2 are 
then (1 : 1) and (2 : 0), the unpolarized and the fully polarized states, respectively. Our 
interest will be in situations when Vc becomes comparable to or greater than the cyclotron 
energy. This again corresponds to > 1, where at z/ = 2, can be seen to be given by 

rs = ^, (2) 

and is clearly a measure of the strength of the interaction relative to the cyclotron energy. 
Here, Landau level mixing becomes crucial and may destabilize the above states. Our prin- 
cipal goal will be to determine the phase diagram of the fully polarized and the unpolarized 
states (2 : 0) and (1:1). 

Our work has been motivated by the recent experiments of Eriksson et al. where they 
investigate by inelastic light scattering both the spin and the charge density collective modes 
at z/ = 2 for samples with rs as large as 6, corresponding to densities as low as p = 0.9 x 10^° 
cm~^. They find a qualitative change in the number and the character of collective modes 
at approximately rs ~ 2. This was interpreted in a Landau Fermi liquid approach, where 
the magnetic field was treated as a perturbation on the zero field Fermi liquid. However, the 
Fermi liquid approach, which is suitable at small magnetic field, is not an obviously valid 
starting point for the problem at hand, and other approaches are desirable. A comparison 
between the ground state energies of the unpolarized and the fully polarized Hartree Fock 
state shows that a transition between them takes place at rs ~ 2.1 for Ez = 0, as we will 
see in Sec.0. This raises the question: Is the ground state at large rs fully polarized? If 



true, this would indeed be an interesting example of an interaction driven ferromagnetism. 
If it is indeed fully polarized, is the (2 : 0) state a reasonable starting point for its study? 
Besides being fully polarized, (2 : 0) incorporates the effect of Landau level mixing, albeit 
in a very special manner, through promoting all electrons in 0| to It- 

A reliable treatment of Landau level mixing lies at the crux of the problem. We shall 
incorporate Landau level mixing in a perturbative time-dependent Hartree-Fock scheme, 
i.e., by incorporating vertex corrections through ladder diagrams in the random phase ap- 
proximation (RPA). The most crucial approximation in our calculations will be a restriction 
to the subspace of a single-particle hole pair; within this subspace, however, the Landau 
level mixing is treated accurately. Clearly, this approach is not quantitatively valid except 
at small r^, but we believe that it gives an insight into the physics even when is not 
small. For a better quantitative description at large r^, it would be important to deal with 
screening by more than one particle-hole pair, but that is outside the scope of the present 
paper. 

As mentioned above, there certainly are parameters for which the fully polarized (2 : 0) 
state describes the actual ground state; in particular, it may occur even when is large 
provided that the Zeeman energy is sufficiently strong. Our approach here will be to take it 
as the starting point and investigate the regime of its stability by calculating the dispersion 
of the charge and the spin density collective modes. The collective modes have a simple 
interpretation when the Landau level mixing is negligible {Vc/hujc 0). The lowest charge 
density excitation mode corresponds to the excitation of one electron from If to 21 Landau 
level. The kinetic energy change of any collective mode is well defined when the Landau 
level mixing is weak, and will be used to label the various modes (this notation will be 
used even when the Landau level mixing is significant, by looking at the evolution of the 
modes from small to large rg). The 11^21 mode will be referred to as the m = 1 mode. 
For spin-density excitation, there are three modes of primary interest, corresponding to 
excitations Oj— >0|, It— and lt^O| which are depicted schematically in (a), (b), and 
(c) in Fig.|l] respectively. The first two are m = modes and the last one is m = — 1 
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mode. These are of course coupled, and a diagonalization of the problem produces the usual 
spin-wave excitation mode, with the energy approaching the Zeeman splitting in the long 
wave length limit, in accordance with the Goldstone theorem, as well as two massive spin- 
density modes. Our principal result is that while the charge density collective mode shows no 
instability, the spin-density collective mode develops a deep roton minimum in the presence 
of substantial Landau level mixing and becomes soft in certain parameter regimes. (Indeed, 
as r 5 is increased, there will eventually also be an instability in the charge density channel, 
indicating the formation of a Wigner crystal, but we have not explored this question since our 
perturbation theory is not valid at very large r^.) Both the existence of the roton minimum 
and its softening as the Zeeman energy is reduced are experimentally testable predictions 
of our theory. The phase diagram thus obtained later is shown in Fig.^]. It is instructive to 
compare it with the phase diagram in the absence of interactions, which consists of the fully 
polarized state (2 : 0) for Ez > Tiojc and the unpolarized state (1 : 1) for Ez < Tiojc with 
a transition taking place at precisely Ez = hue- At small rs, interactions make the (2 : 0) 
state more stable, indicated by the fact that it survives even for Ez < huoc- However, at 
large r^, Ez > huc is required to stabilize the fully polarized state. It is noteworthy that at 
small Ez, the fully polarized state is found to be unstable at arbitrary rs- We also consider 
the collective modes of (1:1), expected to be valid at small Ez and rs- It becomes unstable 
as rs is increased, consistent with an earlier conclusion of MacDonald 0. 

What about the state at large rs but small Zeeman energy? As discovered in our study, 
here the u = 2 state is not described by either of the two aforementioned Hartree Fock 
states. The finite wave vector spin-wave instability of the fully polarized state suggests that 
it is some kind of spin density wave state. It is not possible to be more definitive about it 
based on our present study. Of course, at extremely large rs, when Vc is much larger than 
the cyclotron energy, our calculation is unreliable, but indicates that even if a fully polarized 
state occurs, as expected based on the zero field result, it will most likely not be described 
by the (2 : 0) Hartree Fock state. 

There have been many theoretical studies of the collective excitations of integer quan- 



turn Hall effect (IQHE) states in the past, but, to our knowledge, the spin-density wave 
excitations of the fully polarized (2 : 0) state have not been considered previously. For 
other collective modes, our results reduce to the earlier results in appropriate limits. If 
possible transitions of the electron and hole between different Landau levels are ignored 
and if self-energy corrections are omitted, the problem of collective excitation is reduced 
to determining the binding energy of two oppositely charged particles strictly confined in 
their respective Landau levels f^. In this case the wave function for the bound state is 
independent of interaction potential, and is uniquely determined by the wave vector k. The 
transition of the electron and hole or the recombination of the particle-hole pair has been 
considered in the random phase approximation(RPA) [^0. Later the RPA was incorpo- 
rated with the self-energy correction and the binding energy term by Kallin and Halperin ^ 
where a number of interesting collective excitations from the unpolarized and the partially 
polarized ground state were considered in the absence of Landau level mixing (valid when 
hujc » Vc). The Landau level mixing was considered by MacDonald in Ref 0, treating 
the mixing matrix elements between various modes as small parameters and applying a 
second order perturbation theory. Our calculation will be formulated in terms of diagrams, 
following Kallin and Halperin in Ref. 0, and will be performed with a full treatment of 
the mixing matrix elements within the subspace of a single particle-hole pair at any given 
instant. The diagrammatic formulation of the problem is presented in detail in Sec.|| below. 



In Sec.|T| we describe the diagrammatic formalism used to compute dispersion curves of 
the collective excitation from the fully polarized ground state at z/ = 2. We will concentrate 
on the spin density excitation which will be responsible for an instability of the IQHE state 
in the parameter regime under consideration. Similarly, Sec.^ is devoted to the collective 
excitations of the unpolarized ground state at = 2. Dispersion curves of the spin density 
excitation is computed for various values of rs and Zeeman splitting energy, Ez- The phase 
diagram as a function of and Ez is obtained in Sec.|V] by determining the critical Ez at 
which the energy of the spin density excitation vanishes. From the phase diagram we will 
learn that for large and small Ez neither the fully polarized nor the unpolarized state is 



stable against a spin-density wave state. The paper is concluded in Sec.|V5, where we also 
discuss the implications of our study for the fractional quantum Hall effect (FQHE). 



II. DIAGRAMMATIC FORMALISM OF COLLECTIVE EXCITATION 
A. Algebra in the symmetric gauge 

Even though the choice of gauge does not affect physical quantities, we find it convenient 
to use the symmetric gauge. We will start by establishing the basic algebra in the symmetric 
gauge closely following Ref. [§. 

The Hamiltonian for an electron moving in a two dimensional space under a perpendicular 
magnetic field is given by 

H= — 3 
2m ^ ' 

where the kinetic momentum is written down as 

eA 

7f = -ifiW + — . (4) 

When the magnetic field is uniform, we can use a convenient algebraic method analogous 
to the solution by the ladder operator of the one-dimensional harmonic oscillator. In order 
to construct the algebraic formalism we first note that the x and y components of the kinetic 
momentum are canonically conjugate coordinates: 

T^y] = ■ (V X 1) = (5) 

C Lq 

From the commutation relationship between and Tty we can define a ladder operator so 
that the ladder operator and its Hermitian conjugate satisfy the same commutation relation 
as those of the one-dimensional harmonic oscillator. That is, 

[a,at] = l, (6) 

where 



= ^-^{^x + i7^y)- (7) 
The Hamiltonian can now be written in the form of a one-dimensional harmonic oscillator. 

H = ^{a^a + aa^) (8) 

Therefore the eigenvalues are hujc{n + 1/2) where n is a non-negative integer which is 
known as the Landau level index. The eigenstates, however, cannot be fully determined by 
the Landau level index alone because the energy does not depend on the coordinates of the 
cyclotron orbit center, indicating a degeneracy of the Landau level. Let us define 

(7 = f + ^ (9) 
mujc 

which is convetionally known as the guiding-center operator. The x and y components of the 
guiding-center operator are canonically conjugate coordinates similar to those of the kinetic 
momentum. 

[C,,Cy]^il0' (10) 

Therefore we can define another ladder operator by 

& = + (11) 



which satisfies 



M^l (12) 



and 



[a, 6] = [a^ 6] = [//, 6] = 0. (13) 

The fact that the ladder operator b commutes with the Hamiltonian shows that the degen- 
eracy of the lowest Landau level is actually related to the positioning of the guiding-center 
coordinate. Since we identify two independent sets of ladder operators in a two-dimensional 



|n,m) = i^ilEL|0,0) (14) 



space, the full set of eigenstates can be generated by repeatedly applying raising operators 
on the ground state: 

i\n /, +\m, 

For a magnetic field pointing in the positive z direction, the vector potential in the 
symmetric gauge is given by 

A = ^{-y,x,0) (15) 
and the ladder operators may be written in terms of z{= x + iy) as follows. 

with the Hamiltonian given by 

Here and in the rest of the paper, we use the convention that h = c = e = lQ = l,a.s well as 
the Area= 1. In particular, this implies that the total number of flux quanta piercing the 
system, N^ = ^. 

Up to this point we have concentrated on the single particle Hamiltonian. In order to 
compute concrete physical quantities with the interaction treated perturbatively, we will 
need to use certain matrix elements and their various properties, which we now list 0. 

• Plane Wave Matrix Elements 

The matrix element of a plane wave exp{—ik ■ r) is given by 



where k = y k^'^ + ky^, k = k^ + iky and 

^ ^ 

9npnM = {np\exp{-^Kb^) exp{-^Rb)\na) 

= (|S^)''V-)"^-"'^^:r""(fcV2) (22) 
where, for n^j > Ua, Lni "° is the associated Laguerre polynomial, defined as: 

^"(.)^|:^(r-:) (-) 

For Ufs < Ua define: 

LZ-'^'i^) = ^{-xT'^-^^L^^;-^^ (24) 

Eq.(^T]) can be evaluated first by rewriting k ■ f = [nz + R,z)/2, expressing z and z in 
terms of the ladder operators, and moving all annihilation operators to the right using 

• Matrix Products 

One of the most important properties of the matrix gmm'{i^) is its product 



T.9nA^^l)9ln^i^^2) = 6 '"^'^'Z V.n^ («1 + ^2) (25) 



-K1K2/2 
I 

This can be derived by using the definition in Eq. (^21) and the completeness of the ladder 
operator eigenstates. 

• Eigenfunctions in the Symmetric Gauge 

Similar to duality of the position-space and momentum-space in the one-dimensional 
harmonic oscillator, we have the orbital wave function closely related to the plane-wave 
matrix element. 

{r\n,m) = 0„,m(f) = (-i)"5'mn(^^) /— (26) 

V 2tc 

To derive this, first note that where |0, 0) is given by 

in 



(r10,0) = ^exp[-i^^] (27) 
y 111 4 

which is annihilated by both a and h. The eigenfunction for a general n and m is therefore 
given by 

Now write e~^^/^ = e^^/^e~^^/^ and use 

exp|-l..](| - 2|)-(| - 2|)"exp|i..| = (-2)-''(|)"( A)™ (29) 
Defining t = zz/2 = one gets 

(f|n,m) = — =LJ^=e*/22™z"-'"(— )"t™e-* (30) 



which reduces to Eq. ( p6D with the standard definition of the associated Laguerre polynomial. 

• Trace 

The trace of gmm'{i^) is a constant because the charge density of a completely filled 
Landau level is a constant. 

^ E = ^5fc,o = (31) 

where A stands for the area of system and therefore N^/A is l/27r. 

Using the definition of gmm'{i^), Eq. (l2^) , and properties of the Laguerre polynomial we 
find the transpose and complex conjugate as follows. 

• Transpose 

Quanpi,!^) = gUfjUaij^)- (32) 

• Complex Conjugate 

9nanfi{fi) = gnan^{ — f^) (33) 

• Orthogonality 

The final property we will use is the orthogonality between the plane-wave matrix ele- 
ments which will be useful in computing the self-energy: 

J (fke'''^''^gn^nfX'^)gn^'np'{^<) = 27r5„^,„^/(5„^,„^/ . (34) 
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B. Response function 



The response function is a very important quantity from which we can deduce many 
physical observables assuming hnear response. In particular, the collective excitations cor- 
respond to the poles of the response function. In order to calculate the response function we 
use the zero temperature limit of the Matsubara formalism [|T^ , using the standard analytic 



continuation {iu —>■ u + i6) in order to get the retarded response function. We will describe 
below in detail only the charge-density response function, since the spin-density response 
function can be obtained with a straightforward modification. The charge-density response 
function is related to the density- density correlation function as follows: 

X{k,uj) = - dTe''^^{T^p{k,T)p{-k,0)). (35) 
Jo 

where the density operator, p{k,T), is given in the symmetric gauge by 

The actual computation of the density response function is performed in the perturbation 
theory. Assuming at any given time there is a single particle-hole pair in the system, we 
write the response function which is depicted by the Feynman diagram in Fig.^. 

X{k,iuj)= E (n.,m„|e^^^1n^,m^)D„^„^(^u;)r-';7(fc,2^) (37) 
where the vertex function T^^J^^ {k,iuj) satisfies the following equation: 

X D^^.^^,{tu)TZ:::;\kM 



E E 



xD^^,^^,{zu;)T'::::^,{k,zu). (38) 
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where v{r) is a potential and v{k) is its Fourier transform. The Feynman diagrams cor- 
responding to each term on the right-hand side of the Eq. ( pB]) are shown in Fig.^: The 
first term describes just the bare vertex and the second and third terms depict the RPA 
correction (Bubble diagram) and the binding energy term (Ladder diagram), respectively. 
In the Eq. (^) and Eq. (|3^) Dn^'n^'i^^^) is defined as the frequency integral of the product 
of two Green functions: 

— Gn^/ {iu' - iu)Gnp' {iuj') (39) 

and the Green function is given by 

GnM = , — , (40) 

tuj-{n- i2o)uJc - S;^" 

where n is the Landau level index of the electron or hole and /xq is the chemical potential 
which is halfway between the highest occupied Landau level and the lowest unoccupied 
Landau level. We note that the Green function is fully dressed, containing the self energy 
(S) correction, and is also independent of the m quantum number. In general the self energy 
S is a very complicated function of k and u, but it turns out that in this system the self 
energy is real and depends on just the Landau level index but not k and u. We will use this 
result below, postponing the explicit derivation to the next section. 
Plugging the explicit form of Green function into Eq.(p9|) gives 



9{iio - n^') - 9{^io - rig ) 



(41) 



where 



9{x) 



1 if x > 

if X < (42) 
1/2 if X = 



and it is assumed that the self energy does not modify the Fermi level significantly but 
moves the position of the pole in the denominator of Eg. (|^) . 

In order to solve the vertex equation, Eg. (|38|) , we first eliminate the angular momentum 
index by defining a new vertex function r„^„^(fc, iu;): 

T^^r.,{Ktuj)^e-^'" E g^^^,{-K)K:n;ikM (43) 

t 



Eg .(^41) is obtained from Eg . (1431) through the property of the plane- wave matrix product 
Eg ■ (p5l) ■ Noting that the bare vertex is just a plane- wave matrix element defined in the 
Eg . (^TD and utilizing the properties of the plane- wave matrix elements, we succeed in sum- 
ming over all the angular momentum indices to get a new vertex eguation. 

^rianp ( kjiuj) = — e 9nan0{~l^) 



{k,iuj) 



- E 

ria'np' 

where 

mairifi 



Similarly, the density response function can be written in terms of the new vertex func- 
tion, Tn^n^{k,iu). 

X{k,iuj)= E {-'^T"^''^9n^np{-K.)Dn„np{iuj)Tn^np{k,iuj) (47) 

nari/s 

The third term on the right hand side of Eg.(^) contains a formal expression for 
the Coulomb-potential matrix elements between the wave function ''"^) and 

^ua'nis'i^i^^^) which prove to be the eigenfunctions of the two-body Hamiltonian for op- 
positely charged particles confined in the Landau levels {na^np) and {na'.np), respectively 
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1^. Therefore, following Kallin and Halperin, we will call it the binding energy. The second 
term in Eq.(^) has the physical interpretation of the exchange energy of particle-hole pair 
which is the only term in the usual RPA. For convenience let us define riap = — np, and 
denote the binding energy and the exchange energy as Vnau"/ (k) and Un^nf (k) respectively: 

K:n/ik) ^ r"''-""^'e-'=^/^^„„„,(-/.)^^„,,., (-«:) (48) 
VZi;/{k) ^2nJ d'fl J d'm:^^^^{fl,rMri - r^Wn^, ,nAr-{,T,) (49) 



Then Eq. becomes 



2"°" _l..2, 

E - U:i-/{k)]Dr,^,^^,{iu)T^^,^^,{kM- (50) 



The Green function Dn^^n/j contains a factor [OiriQ + 1/2 — rip) — 9(nQ + 1/2 — n^)^ which 
vanishes for certain choices of {rtainp). However, these choices co not contribute to x iii 
Eq. (|47|). Therefore, restricting to only those values of («„, ra^) for which Dn^ n^ is non-zero, 
we transform the equation into an associated set of linear equations: 

= '^e-''"9n.n,{-Kl (51) 

where 

D-\nA'^) = {luJ - {< - V)^c - iKl' - Kid) X {e{n, + 1/2 - V) - ^^(^0 + 1/2 - <)). 

(52) 

(Note that here D~\^^,{iuj) is the inverse of the matrix element, not the matrix element of 
the inverse.) Eq . ([5l|) can be written in the form of a matrix equation by defining a matrix 
M as follows: 

M::if{k,tu;) ^ 5„.,n„'5„„„,'/^;X.(^a;) + VZi^^' {k) - U^'Z/ik). (53) 



Mn2np {k,ico) can be viewed as a matrix element of M if we consider the set of indices 
{iT-afifs) to be a collective index. Therefore we write Eq.(|^) in the following form: 

Y: MZ'::f{kMDn^'nA^uj)T^^.^^,{kM = -_e-^"/2-^^„^ (-«:). (54) 

na'rif}' 

Inverting the matrix M amounts to solving the vertex equation. That is, 

where is the inverse matrix of M. Finally, we obtain the density response function in 
terms of the matrix M. 

x{kM=i: E e-'^'/'gr,^^^{-R){M-X::;ikM9n.'nA-'') (56) 

nani3 na'np' 

Since the summation of the Landau level indices should be performed over all the possible 
states, the number of terms is actually infinite. However, it is possible to obtain an accurate 
estimate keeping a reasonably small number of terms. 



C. Collective excitation 

The collective excitations are the poles of the response function, which, from Eq.(^ 
correspond to energies for which the inverse of M becomes singular, that is to say, when 



Det 



M{iuo ^uo + i5) =0. (57) 



We carry out detailed computations of the binding energy, RPA energy and self energy 
in order to explicitly evaluate the dispersion curve of the collective excitation. We start with 
the binding energy. 

• Binding energy (Ladder diagram contribution) 

The explicit form of ^naU/ji'^i-i^^) is obtained from Eq.(^) by utilizing the plane-wave 
matrix product formula, Eq. fpSl) : 



27T 



mamp 



-i)"'"'^ 2 2 2 



27r 



-I) 



-k'^/i -{ri'^+r2^)H iKZ2/2 



27r 



-A:2/4^-(ri2+r22)/4j(K22+K2i)/2^fi 2:2/2 



2tt 

■ \n. 



e^'^"'^ 9n^ngi-K. - izi + iZ2) 



JR-{k+rxz/2) 
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_ |r+fcxg| 

e 4 



9n^n0{-i{z + ik)). 



(58) 



Following Kallin and Halperin, the wave function ^'JJ^ "^2) can be shown to be an 
eigenstate of the following Hamiltonian where the particles 1 and 2 are projected onto the 
Landau levels with indices na and respectively. 



H = P HP 



Pr 



a, '1-13 



-> 2 -> 2 

2m 2m 



Pr 



hUciria + 12/3 + 1) + Pn^,npVarb{\ri - r^|)-Pn,,n^ 



(59) 



where vf^ (vf^) is the kinetic momentum for the particle 1 (2), farfe('") is an arbitrary (at- 
tractive) potential and Pna,np is the operator projecting the particles onto the Landau levels 
with indices ria and n^. One can prove that ^^(^'l, t^) is the eigenfunction of H by 
computing the overlap between farfe(|n — ?^|) x *&n„,n^(n, ''''2) and any arbitrary basis state 
in the projected Hilbert space, for example (j)nam^{ri)4>n0ms{f^) ■ Then we note that it is 
proportional to the overlap between $5^^^„^(?''i, "^2) and 4>ncm^{f\)<Pnf,ms{^)- Of course, the 
proportionality constant is the eigenvalue which is equal to the binding energy previously 
defined in Eg. (^9]) in the case of Ua = and = n^. 

Now let us get the explicit formula for the binding energy. 

V^ifik) = 2n f d\-l f rfV2$:;^,„^(rI,r-2)t;(rl -r-2)$:;^,„^,(rl,r-2) 



■afJ — "-al3 



2ti 



d R d re 2 y(^r)g^ (-{(^z + ik))gn^>na'{,-i{,z + ik)) 
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y2"«2"" na\na \ ) J (27r) 

where is an associated Laguerre polynomiaL By using the fact Jq^ d9e^^^^°^^~^'^^^ = 
2ni"'Jn{x) we perform the angle integrations to get the final formula: 

2""2"" ra„!n„'! / io 27r ^ 



X I rfre-^/V"'^''+"'^''+^L::;''(-)L:y (-)J„^,_.^,,(gr), (61) 

where no/3 = na — Ujs, Uap = nj — np and J„ is a Bessel function and is the angle of k 

measured from the x-axis. 

• Random phase approximation energy (Bubble diagram contribution) 

The exchange energy from the RPA is rather straightforward to calculate since there is 

no integration involved. The explicit form is given by 

Ki::\k) = ^""''-""'''e-'^'/^^7n.„,(-«:)^^7„.'„^(-«:) 

l2n.2«<^'njn,'!j ^ ^n, (y)V ^^^^ 



Incidentally, a comparison between Eg. (|6T1) and Eq.(|62D reveals that Un^nf {k) and 
Vn^nl'^ (k) have the same phase factor, e*^""''"""''')^^. Therefore the phase factor can be elim- 
inated in a consistent way, which is expected because the system is uniform and isotropic. 

• Self energy 

The last diagram in the Fig.^ is the self-energy contribution to the "full" Green function. 
As before, a single particle-hole pair is assumed, which is reflected in the Feynman diagram 
through the unscreened Coulomb line. The corresponding Dyson equation is solved in the 
conventional way: 



1 

iuj - {ua - /io)t^c - S;^o 
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(63) 



where 



G^^\iu) = , — 

«cu - K - /io)cuc 



(64) 



and 
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In the above equation the exphcit form of the single particle eigenstate, Eq.(^), and the 
plane- wave matrix product formula, Eq . (|25|) , have been used. It is convenient for the com- 
putation of integrals to change the variables from f\ and to the center of mass coordinate 
R and the relative coordinate r. In the form of a complex number R ^ Z = {zi + Z2) and 
r z = Zi — Z2. 
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The angular momentum index ma is naturally eliminated after the integration over the 
center-of-mass coordinate by using Eg. (l3^) . The plane- wave matrix product formula has 
also been used when we proceed from the third step to the fourth step. Using J2n=o -^n(^) = 
-^no '^^^ finally write down the self energy as follows: 



In order to see what modifications need to be made in the above analysis to include the 
spin degree of freedom, it is instructive to recall the physical meaning of the three parts of 
the collective exciton energy. The binding energy due to the ladder diagram is the direct 
interaction between the excited electron and the hole. Therefore it will not be affected by the 
presence of the spin degree of freedom. On the other hand, the RPA energy is the exchange 
energy between the excited electron and the hole. So it will vanish if the excited electron 
has a different spin than the hole, as in the case of the spin density excitation. In other 
words the electron-hole pair with the same spin cannot recombine through the Coulomb 
potential. In the case of the charge density excitation, however, the RPA energy depends 
on the polarization of the ground state. We have computed the RPA energy in the previous 
section assuming that all the electrons have the same spin, which corresponds to the fully 
polarized state. If the ground state is unpolarized, the RPA energy will be twice as large 
as that of the fully polarized state simply because the particle-hole pair can be created and 
annihilated as either a spin-up or spin-down pair. Formally speaking, the vertex equation. 




(67) 



D. Spin degree of freedom 



Eg. (pHD , will have two identical RPA terms for a given set of indices. The self energy term is 
due to the exchange energy between a given electron and the rest of electrons in the system 
while there is no direct term because we assume a neutralizing positive charge background. 
Since it is an exchange term, we have to include only the interaction between the electrons 
with the same spin. Therefore Eq.(|67|) is generalized to include the spin degree of freedom 
as follows: 



2 ^2 



/"OO ft roo f-^ fp-^ ^ 

^l^''^ = -l dr-rLl{-)L\^^{-)Uqr)e-^''\ (68) 



where 



1/2 for the spin antiparallel to B field 
a={ ' (69) 

-1/2 for the spin parallel to B field 

and n{a) is the index of the highest Landau level occupied by the electron with spin a. 
Using this new self energy one can rewrite Eq.(^T]) as follows: 

n (■ ^_ g(/i(ff^) - np) - e{ii{ao) - n») 

where the chemical potential yu(o") was defined earlier, and the Zeeman coupling is included 
through the term aE^. 



E. Solutions for the pole of the response function 

Computing the dispersion curve of collective excitations for a general rs requires solving 
the equation for the pole of the response function, Eq. (|57|) . In the limit rs 0, when there 
is no Landau level mixing, the matrix has a block diagonal form since collective modes with 
different kinetic energies do not couple, and each block, which has a finite dimension, can 
be diagonalized separately to obtain the collective mode energies At non-zero values 
of rs, however, the full matrix must be diagonalized. Strictly speaking, the matrix M in 
Eq.(^7l) is of infinite dimension, but in practice, we work with a finite size matrix, keeping 
a sufficient number of Landau levels to ensure a convergence of the collective mode energy. 
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For the lowest energy collective modes, which are our primary concern, and for < 6, we 
find that it is adequate to work with M of dimension of up to 20. 

We also find it useful to convert the Eq. (0) into an eigenvalue equation, the solution 
of which can be obtained using standard linear algebraic methods. In this formulation it is 
natural to define an effective Hamiltonian to be the matrix in the eigenvalue problem. The 
details of the procedure are discussed in the remainder of this section. 

For convenience we write down Eq. (|56|) here after the analytic continuation. 

x{kM=Y. E ^-^^ e-'=^/Vn,(-/^)(M-^)::l7'(A:,c^)Wn,'(-'^)- (71) 

For an arbitrary set of indices (n^, np) where is greater than np, there is a reversed set 
(ripjUa) whose kinetic energy cost is negative. The mode described by the reversed set of 
Landau level indices was called a negative-energy mode by MacDonald because of its negative 
kinetic energy cost. MacDonald considered mixing between positive and negative energy 
modes as well as between the positive energy modes in second-order perturbation theory in 
order to compute the collective excitations at general rs in the spin unpolarized ground state. 
In the present work we directly solve the pole equation, Eq. f^Tf) , instead of approximating 
it to the second-order. In computing the lowest lying mode it is especially important to 
consider mixing with the negative-energy modes since the lowest positive-energy mode is 
energetically closest to the negative-energy modes. 

In order to make our discussion concrete and transparent, let us explicitly write down 
the matrix elements in the case of the full spin polarization. If — np = m > 0, n^' > /ig 
and np < /iq, Eq.(0) is written, after the analytic continuation, as follows: 

Kl'uA^) = - muc - iKl' - Kl'))- (72) 
If the order of the indices is reversed, we get 

D-}nA^) = {uJ + muc + {Kl' - Kl')) X (-1) 

= D~\^^,{-u) (73) 
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Therefore the matrix M defined in Eq.(Q) and its counterpart with reversed indices are 
and 

where we set = without loss of generahty and therefore V^I^n^" {Un^rZ" ) = 

[ I) 1^ Vnari/s {Urianfj )■ 

The sign in front of the second term in Eq.(^) can be interpreted so that the negative- 
energy mode has the angle of wave vector equal to vr: for the negative-energy mode 6'^ = vr in 
the phase factor of Kl)fn"° and Unj^rZ" ■ According to Eq. ([75|) the negative-energy mode has a 
negative frequency and the opposite direction of wave vector relative to the positive-energy 
mode. Therefore when the positive-energy mode is viewed as a plane wave e**^"^*"^'^, the 
negative-energy mode is the complex conjugate plane wave e~*(^*~'^'^. In this interpretation 
the requirement of negative-energy mode is natural since an arbitrary plane wave with 
k is written as a linear combination of e*'^'^*"'^'^ and e~*('^*~'^''^. Therefore in general a 
collective excitation with k should be described by both the positive-energy and negative- 
energy modes. Incidentally, we mention that the mode describing an excitation within the 
same Landau level, for example the spin-wave Goldstone mode, does not have the negative- 
energy counterpart because there should not be any double counting in Eq. (|7TD . 

In any case we reahze that the pole equation, Eq.(^), is not an eigenvalue equation 
as it stands because of the sign change in u for the negative-energy mode. But it can be 
transformed to an eigenvalue equation as follows. Let us denote M in terms of sub-matrices. 



M(fc,cj) 



(76) 



Moo(fc,c^) Moi(/c) 
Mio(A;) Mn(A;,u;) 

where Moo(Mii) is associated with mixing between the positive-energy (negative-energy) 
modes and Moi(= Mio) is between the positive and negative energy modes. Thanks to 

9?? 



M{k,uj) 



(78) 



a property of determinant we can obtain the solution of Det[M{k,uj)] = by solving the 
following equation. 

Det[M{k,uj)] = 0, (77) 

where 

Moo(fc,cu) Moi(/c) 
-Mw{k) -Mii{k,uj) 
Now one can define an effective Hamiltonian matrix H{k) using M{k,uj). 

M{k,uj) = H{k) -ul (79) 

Therefore Eq.(|77|) amounts to the eigenvalue equation of the effective Hamiltonian matrix 
H, which, however, is not a Hermitian matrix because of the sign change. 

Solving the eigenvalue equation of a non- Hermitian matrix is complicated by the fact that 
the eigenvalues of a non-Hermitian matrix can be highly sensitive to small changes in the 
matrix elements . The sensitivity of eigenvalues to rounding errors during the execution 



of some algorithms can be reduced by the procedure of balancing. The idea of balancing 
is to use similarity transformations to make corresponding rows and columns of a matrix 
have comparable norms, thus reducing the overall norm of the matrix while leaving the 
eigenvalues unchanged. Then the general strategy for finding the eigenvalues of a matrix is 
to reduce the matrix to a simpler form, and perform an iterative procedure on the simplified 
matrix. The simpler structure we use is called the Hessenberg form. An upper Hessenberg 
matrix has zeros everywhere below the diagonal except for the first subdiagonal row. Then 
one can find the eigenvalues by applying the QR algorithm repeatedly to the Hessenberg 
form until convergence is reached. 

Finally, we mention that even though the above formalism has been developed for a 
general situation, we will concentrate on the physics at z/ = 2 in the following sections 
taking the Coulomb potential as the interaction. When the finite thickness effect of the 2D 
system is of interest, one can replace the Coulomb potential by an effective potential such 
as the Stern- Howard potential [|l^,|ll|]. 
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III. COLLECTIVE EXCITATIONS OF THE FULLY POLARIZED IQHE STATE 

AT v = 2 



Equipped with the exphcit formulas for the binding energy, the RPA energy and the 
self energy, we compute the dispersion curves of the collective excitation from the fully 
polarized ground state. First we study the large B field limit, i.e. the small rs limit, 
where the (time-dependent) Hatree-Fock approximation is valid. In the small rs limit the 
effective Hamiltonian H defined in Eq. ([79|) is already block-diagonalized so that only the 
matrix elements within the Hilbert subspace of the same kinetic energy survive. The off- 
diagonal terms due to the interaction energy become negligible compared to the kinetic 
energy. Therefore the pole equation is simple to solve in this case. For an arbitrary value of 
rs the Hamiltonian is generally complicated and needs to be diagonalized over whole Hilbert 
space. 

When we consider the small rs limit, the integrals we encounter in Eq. (|6ll) , Eq.(|6^) and 
Eq.(|67|) can be expressed in terms of /„(«,/?) and defined as follows: 



where n is an integer. The explicit functional forms of fn{a,P) and (7„(a;) are categorized 
into ones with even n or with odd n: 




(80) 



and 




(81) 




(82) 



where 




/i(«,/5) 



2a 



(83) 



and Iq is a modified Bessel function. Similarly, 



92m+l(0l) = 



da ^ 
da 



(84) 



where 
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A. Charge density excitations 

Since the lowest mode in the energy spectrum is most relevant, we consider the mode 
where an electron is taken from the {n — 1) Landau level and promoted to the {n — 2) 
Landau level without flipping the spin in order to obtain a charge density excitation. Since 
its kinetic energy is huj, we will call the mode the m — 1 mode. As discussed before, this 
mode is not mixed with other modes in the limit of small rs, i.e., the matrix M is already 
diagonal. For the lowest charge density excitation we have the following equation: 

Ml'^,{k,uj) =uj-ujc- ^(Si„=2 - - ^Ulj{k) + ^V^:hk) = 0. (86) 

In the above equation we defined dimensionless matrix elements so that 

= cir'(^)/(^v^^o) 

Kl = Kl/{e'/el,) (87) 

Let us denote the solution of the pole equation as A(A;) from now on. Then the dispersion 
curve of the lowest charge density excitation is given by 

9fi 



e2 



^^c + -r 



92{a) - ^^4(«) + T^geia) 
Z Lb 



a=l/2,l3=k 



.£e--.(l-g (88) 

As explained earlier, Landau level mixing in the non-zero rs regime is included by di- 
agonalizing the effective Hamiltonian defined in Eg . (|79D . Using the term m = 1 mode for 
the sake of convenience to indicate the lowest charge density excitation in the fully polar- 
ized state at a general rs, we plot their dispersion curves in Fig.^ which shows that the 
charge density excitation modes do not exhibit any sign of an instability in the parameter 
range of rs considered here. A non-trivial check of our calculations is to make sure that the 
collective excitations computed within our approximation satisfy the exact Kohn's theorem 



which states that the m = 1 mode energy must approach tuc as /c — Fig.^ shows 
that Kohn's theorem is satisfied not only for the pure mode but also the mode with Landau 
level mixing in general rs- Incidentally, the Zeeman coupling will not affect the dispersion 
curves because the spin is not flipped. 



B. Spin density excitations 

Following the convention used in the previous sections, we indicate the excitation modes 
in terms of m, the kinetic energy of the mode in units of huc in the limit of small rs- We 
shall see that, unlike the charge density excitation, the lowest lying spin excitation can be 
either m = —1 mode or the lower one of the two m = modes, depending on the value 
of rs- The m = — 1 mode describes the process whereby an electron in the n = 1 Landau 
level is demoted to the n = Landau level with its spin reversed, whereas the m = mode 
has an electron with its spin flipped in the same Landau level. In the latter case, there are 
two possible modes: spin density excitation within n = 1 Landau level or n = Landau 
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level. At small the m = —1 mode the lowest excitation while for large the lowest spin 
density excitation is a m = mode. Also, without the Zeeman coupling, the spin reversed 
mode always causes an instability of the fully polarized (2:0) state. A determination of the 
Zeeman splitting energy Ez required to make the (2:0) state stable for general rs will be 
one of the main goals when we try to obtain the phase diagram of the spin polarization as a 
function of and Ez- The dispersion curves of the pure spin density excitation, however, 
can be computed in the limit of small rs without recourse to the actual value of Ez- We 
assume that it is large enough to stablize the excitation because the Zeeman energy is just 
a constant shift in this limit. 
• m = — 1 mode 

As in the case of the charge density excitation, we first solve the pole equation in the 
small rs limit to get the pure mode without any Landau level mixing. With matrix elements 

^2 



Mqi (/c, cu) = uj + ujc 







the solution is 



e/o L 
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and 



-V,y{k)=[~foia,P) + -Ma,/3) 



2V2 



a=l/2,/3=k 



;i + y)/o(^)-y/i(x: 



(92) 



The dispersion curve of this m = — 1 mode is plotted in Fig.^ 
• m = modes and the spin wave mode 

Since there are two possible m = modes, the pole equation becomes a matrix equation 
even in the small rs limit: 
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The matrix elements are given by 
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and 
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Since spontaneous symmetry breaking occurs in the fully polarized state, there must 
be a spin-wave Goldstone mode whose energy approaches the unshifted Zeeman splitting 
energy in the long wavelength limit. Similar to the case of the charge density excitation 
it is important to check if the spin density excitations computed within our approximation 
satisfy this exact theorem. When we take the k ^ limit of Eg . (|9^) , we have the following 
equation: 
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The solutions are 



Ai(A; = 0) = (98) 

^^(*= = °' = £/|- 

This confirms the existence of a Goldstone (spin-wave) mode. Furthermore it can be shown 
that the massless spin-wave mode exists with an arbitrary potential and Landau level mixing 
within our approximation. The dispersion curves of the two m = modes are also plotted 
in Fig.^ along with that of m = — 1 mode. 

As before, the dispersion curves for general and Ez are obtained from diagonalization 
of the effective Hamiltonian. We take the dispersion curves at rs = 1.0 and 3.0 as examples 
and plot them in Fig.^ to illustrate the qualitative difference between the small and large 
rs regimes. We will sometimes use the term Goldstone mode in order to indicate spin-wave 
excitation mode for general rs since its energy approaches Ez as k goes to zero. We will 
also use the term m = —1 mode for the lowest excitation for small rs since it has the lowest 
energy in the limit of vanishing rs- The physical implications and the corresponding phase 
diagram will be discussed in more details in the later section. 



IV. COLLECTIVE EXCITATIONS FROM THE UNPOLARIZED IQHE STATE 

AT u = 2 

The collective excitations in the spin unpolarized state are computed in this section in 
complete analogy with the previous section. First, let us consider the small rs limit with 
zero Ez- The energy of the lowest charge density excitation is 



do elo 



(100) 



a=l/2,f3=k 



And the energy of the lowest spin density exciation is 



A(fc)=o;c + ^P°=i-SLo] 



el 



elo 



(101) 



J o=l/2,/3=fc 

The dispersion curves of the above pure modes are plotted in Fig.^. We can use Kohn's 
theorem to check if our approximation is reasonable. Fig.|^ shows that the energy of charge 
density excitation approach tuc as A; — > 0. Since the charge density excitation does not show 
any instability, we will concentrate only on the spin density excitation to obtain the phase 
boundary which is the critical Zeeman splitting energy needed for the stable excitation. 
Fig.|| shows the dispersion curves of the lowest spin density excitation for r5 = 1.0 and 
various values of Ez- 



V. THE PHASE DIAGRAM 

The phase diagram of the state at z/ = 2 can be obtained at two levels of sophistication. 



A. Phase diagram in Hartree Fock approximation 

The simplest approximation is that of non-interacting electrons. In this case, the phase 
boundary is given simply by Ez = huJc, as shown in Fig.^. As we shall see, this is sensible 
only in the limit of rs 0; interactions modify the phase diagram substantially elsewhere. 

In the simplest approximation, interaction can be incorporated by comparing the energies 
of the fully polarized and unpolarized ground states in the Hatree-Fock approximation, that 
is, by assuming that the ground state contains either Of and OJ, Landau levels fully occupied 
or Ot and If. The contributions of the kinetic energy and the Zeeman coupling to the 
ground state energy are straightforward. The exchange interaction energy in the Hartree- 
Fock approximation can be evaluated in terms of the self-energy defined in the previous 
section. The self-energy S^*^"") is the exchange interaction between an electron in the Landau 
level with index n and all the other electrons with the same spin. The exchange interaction 



energy per particle is then the sum of the self-energies for all electrons divided by two times 
the number of electrons where the factor of two prevents a double counting. That is to say, 



1 



1/2 n(a) 



^- = ^ E (Esr)- (102) 

a=-l/2 n'=0 

Therefore the ground state energy per particle is 

E, = l E [K + ^)^^c + aE^ + ^S:^)]. (103) 

a=-l/2n'=0 ^ ^ 

The average energy of the fully polarized ground state is 

Similarly the average energy of the unpolarized ground state is 

E„(l : 1) 1 1 ~n 1 /tT , , 

The phase boundary is given by the solution of the following equation: 

E,{2 : 0) - E,{1 : 1) _ 

hue 

Therefore the critical Ez/huc as a function of rs is 

Ez , 3 /7f 



0. (106) 



shown in Fig.^. 

B. Phase diagram from collective mode instability 

The phase diagram obtained by a comparison of the energies of the Hartree Fock states 
is not fully reliable for two reasons. First, it neglects Landau level mixing, which is crucial 
for the issue of interest here. Secondly, it does not allow for the possibility of other states 
in the phase diagram. Therefore, it is more appropriate to look for instabilities of the two 
states by asking when one of the collective modes becomes soft. Indeed, a first order phase 



transition may (and likely will) occur even before the collective mode energy approaches 
zero, but we believe that the phase diagram obtained by considering instabilities ought to 
be reliable qualitatively and even semi- quantitatively. 

As noted previously, there is no instability in the charge density wave collective mode 
in the parameter range considered here. We have determined the onset of the spin density 
collective mode instability both in the fully polarized and the unpolarized states by varying 
rs and Ez, as shown in Fig.^ and Fig.^]. Fig.^ shows the phase diagram thus obtained. The 
following features are noteworthy. 

(i) The nature of instability is different depending on whether rs is small or large. At 
small rs, where the interactions are negligible and the physics is dictated by the Zeeman 
energy, the lowest energy spin-density excitation is clearly the m = — 1 mode. This continues 
to be the case for ^ 2; here the m = — 1 spin-density mode is responsible for the instability 
of the (2:0) state, as shown in Fig.^. However, for rs ^ 2, interactions are sufficiently strong 
that the spin-wave mode becomes the lowest energy mode and causes the instability. One 
way to understand why the m = spin wave mode has lower energy at large rs than the 
m = — 1 mode is because whereas the former has an energy equal to Ez in the long wave 
length limit no matter what r5, as guaranteed by the Goldstone theorem, the energy of the 
latter is determined by the interactions. 

(ii) At small rs, the interactions make the fully polarized (2:0) state more stable as 
compared to the non-interacting problem, as evidenced by the fact that the transition out 
of it takes place at a Zeeman energy smaller than hujc- This is precisely as expected from 
the exchange physics, as also captured by the Hartree-Fock phase diagram. 

(iii) The instability in the spin-wave mode occurs through the development of a roton 
minimum, the energy of which vanishes at certain Zeeman energy. The roton minimum in 
turn is caused by the Landau level mixing, underscoring the important role of Landau level 
mixing at large rs- Without Landau level mixing the spin- wave mode does not show any 
instability as shown in Fig.^. 

(iv) For the unpolarized state, the lowest spin density excitation is the m = 1 mode at 



all rs, the long wavelength limit of the excitation energy of which is fixed to be hujc ~ Ez 
independent of r^. Here also, the instability occurs through a roton, which becomes deeper 
as the is increased (and interactions become stronger). Since a spin flip is favored by 
exchange, the energy of the m = 1 mode decreases with increasing rs, consistent with the 
feature that the critical Ez in this case is monotonically decreasing as a function of rs, as 
is shown in Fig.|[ 

(v) At small rs there is a small region in Fig.^ where the two phases coexist. This 
clearly is an artifact of various approximations involved in our calculation, and implies that 
the actual locations of the phase boundaries cannot be taken too seriously. As mentioned 
earlier, a first order transition is likely to occur before the roton energy vanishes. 

(vi) At large rs and small Ez, the ground state is derived neither from (2:0) nor from 
(1:1). Since we find a finite wave vector instability in the spin density wave excitation, it 
is natural to expect that the state here has a spin-density wave state. Further work will be 
required to establish the nature of this state in more detail. 



VI. CONCLUSION 

The principal outcome of our calculations is the phase diagram in Fig.^ which shows 
the regions where the fully polarized and the unpolarized Hatree Fock states (2:0) and (1:1) 
are valid. We believe that it should be possible to investigate the roton minimum in the 
spin-wave excitation of the fully polarized state as well as its instability in inelastic light 
scattering experiments. 

Another situation where similar physics may apply is in the case of composite fermions 



T5| at effective filling factor z/* = 2, which corresponds to the electron filling factor v = 2/5. 
It is easier in this case to see a transition between the fully polarized and the unpolarized 
states because the effective cyclotron energy for composite fermions is substantially small 
compared to the cyclotron energy of electrons, which makes it possible to obtain Ez com- 
parable to or larger than the effective cyclotron energy in tilted field experiments. The 



critical Zeeman energy at the transition was calculated by Park and Jain by comparing the 



ground state energies |T^, in reasonable agreement with the experiments of Kukushkin, von 
Klitzing, and Eberl ||18|. Park and Jain also estimated a mass for composite fermions, the 
"polarization mass" by equating the critical Zeeman energy to the effective cyclotron energy. 
There is one subtlety though. In the case of composite fermions both the effective cyclotron 
energy and the effective interactions derive from the same underlying energy, namely the 
Coulomb interaction between the electrons, and therefore neither the interactions between 
composite fermions nor a mixing between composite fermion Landau levels can, in principle, 
be neglected in any realistic limit. These would provide a correction to the mass obtained 



in Ref. |T^. However, we note that the mass was reliable to no more than 20-30%, and the 
corrections may be negligible compared to that. 



Interestingly, there is experimental evidence [p!7| , p!8| that the transition between the fully 
polarized and the unpolarized composite fermion states does not occur directly but through 
an intermediate state with a partial spin polarization. Murthy |]19[ has proposed that this 
state is a Hofstadter lattice of composite fermions, and has half the maximum possible 
polarization. It would be interesting to see if similar physics obtains for z/ = 2 as well. In 
particular, the phase diagram of Fig.^ predicts that for > 1, the transition from the fully 
polarized state to the unpolarized state as a function of the Zeeman energy is not direct but 
through another, not yet fully identified state. (We suspect that this may be true at any 
arbitrary rs, although not captured by our calculated phase diagram.) 

This work was supported in part by the National Science Foundation under grant no. 
DMR-9986806. We thank G. Murthy for discussions. 
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FIGURES 
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FIG. 1. Schematic diagrams for some relevant excitations contributing to the lowest spin density 
excitations from the fully polarized IQHE state at = 2. The Zeeman splitting is set to zero for 
simplicity. The processes of moving an electron from 0| to OJ, Landau level and from 1| to 1| 
Landau level are described in figure (a) and (b) respectively. Figure (c) depicts the process of 
moving an electron from 1| to OJ, Landau level. These are labeled by the kinetic energy change in 
the small rs limit as m = (a and b) and m = — 1 (c) modes. 




FIG. 2. Phase diagram of the u = 2 state as a function of Ez/ijoc and rs- The phase 
boundaries are computed from the roton instabihty of the spin wave excitation of the fully polarized 
and unpolarized state. The region where the fully polarized state is stable is denoted by (2:0) while 
the region for the unpolarized state is denoted by (1:1). 




FIG. 3. Feynman diagrams for the response function xik,u}), vertex function Tuan/ {k,uj) and 
the self-energy correction to the "full" Green function G„„ (a;) . 
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FIG. 4. Dispersion curves of the lowest charge density excitation from the fuUy polarized 
IQHE state at = 2 for various values of rs- ^Prom the top, the values of rs are 0.0, 1.0, 2.0 and 
3.0. We denote this mode by m = 1 charge density mode, since its kinetic energy approaches unity 
(in units of huic) in the limit of rs = 0.0 there is no Landau level mixing. 
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FIG. 5. Dispersion curves of the spin density excitations in the fully polarized IQHE state 
at u = 2 for rg = 0. Note that when = the Zeeman energy contribution is a constant shift 
in energy. Each mode of the spin density excitation is denoted by m i.e. its kinetic energy in 
units of hue- Since there are two degenerate modes for m = 0, we diagonalize the Hamiltonian 
in this Hilbert subspace to get the eigenstates whose dispersion curves are plotted as the thick 
solid and dotted lines in the top graph. The thin lines in the top graph indicate the uncoupled 
modes described in the adjacent diagrams where j (j) denotes the spin-up (-down) electron and 
the horizontal lines denote the Landau levels. The lower mode in the two m = modes is referred 
to either as the Goldstone mode or the spin- wave excitation mode. The dispersion curve of m = — 1 
mode is plotted in the bottom graph. 




FIG. 6. Dispersion curves for the spin wave excitation and the m = — 1 spin-density mode in 
the fully polarized IQHE state aX u = 2 for = 1.0 and 3.0. The solid line corresponds to the 
spin wave excitation while the dashed line to the m = — 1 mode. The various dispersion curves 
are plotted for different ratios of Ez/hu)\ in ascending order of line thickness, this ratio is given 
by 1.2, 1.0, 0.8 and 0.7 in the top graph and by 1.8, 1.6, 1.4 and 1.2 in the bottom graph. The 
kinks in the dispersion appear at the anticrossings between different modes. (Only the lowest two 
collective modes are shown here.) While the m = — 1 mode has lower energy at = 1, the spin 
wave mode is the lowest energy mode for = 3. 
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FIG. 7. Dispersion curves of the lowest charge (dashed hne) and spin (sohd hne) density 
excitation in the unpolarized IQHE state at u = 2 aX rs = 0. The Zeeman energy has been set to 
zero in the plot; for rs = the Zeeman energy contribution to the spin density excitation dispersion 
is a constant shift of Ez- The charge density excitation is not affected by the Zeeman coupling. 
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FIG. 8. Dispersion curves of the lowest spin density excitation in the unpolarized IQHE state 
at 1/ = 2 for rs = 1-0 and various Ez^s. In ascending order of hne thickness the values of Ez/tiu 
are 0.50, 0.60, 0.65 and 0.70. 




FIG. 9. Phase diagram of the v = 2 state as a function of Ez/toc and rs- The soUd 
Unes indicate the phase boundaries which are computed from the roton instabihty of the fully 
polarized and unpolarized states. The dashed line indicates the phase boundary obtained from 
the comparison between the two ground state energies in the Hartree-Fock approximation. The 
long-dashed line is the phase boundary of non-interacting electron system. The region where the 
fully polarized state is stable is denoted by (2:0) while the region for the unpolarized state is 
denoted by (1:1). 
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